00 

o 

O 



O 



c^ 



The quantum-limited comb lineshape of 

a mode-locked laser: Fundamental limits 

on frequency uncertainty 



J. K. Wahlstrand,' J. T. Willits,' ^ q r. Menyuk,^ and S. T. Cundiff' 

^JILA, National Institute of Standards and Technology and the University of Colorado, 
Boulder, CO 80309-0440 
-Department of Electrical and Computer Engineering, University of Colorado, Boulder, CO 
CN ■ 80309 

^Computer Science and Electrical Engineering Department, University of Maryland 
Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250 



*vq ' Abstract: We calculate the quantum-limited shape of the comb lines from 

a mode-locked Ti: sapphire laser using experimentally-derived parameters 
^.' • for the linear response of the laser to perturbations. The free-running width 

C_) I of the comb lines is found across the laser spectrum. By modeling the 

effect of a simple feedback loop, we calculate the spectrum of the residual 
M-(| phase noise in terms of the quantum noise and the feedback parameters. 

Finally, we calculate the frequency uncertainty in an optical frequency 
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^ I measurement if the limiting factor is quantum noise in the detection of the 

. ,-H ' optical heterodyne beat. 
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1. Introduction 

When light passes though any gain medium, it will acquire noise through spontaneous emission 
(SE) m. The SE noise sets a fundamental limit on the linewidth of a laser. For a continuous 
wave laser, this Schawlow-Townes limit is due to phase jitter Q. The situation with mode- 
locked lasers is more complicated. A mode-locked laser produces a train of pulses, regularly 
spaced in time, and the frequency spectrum is a series of narrow comb lines. The width of 
each comb line depends on the timing and phase jitter of the pulse train. The laser dynamics 
are characterized by four pulse parameters: the pulse energy, the central or carrier frequency, 
the central pulse time, the phase, and in addition a fifth parameter: the round-trip gain. In all 
mode-locked lasers, strong nonlinearities couple amplitude and frequency fluctuations into tim- 
ing and phase jitter The quantum-limited noise properties of mode-locked lasers were first 
treated by Haus and Mecozzi in |[3l, hereafter referred to as HM. Until recently, subsequent 
work focused mainly on quantum-limited timing jitter With the development of femtosecond 
frequency combs ID, the more general problem of timing and phase jitter has received more 
attention |l5]|6l|7l[8l|9l[l0l[n][Ill- While HM calculated the timing and phase jitter rather than 
the comb linewidths, that paper contains almost all of the information needed to calculate the 
line widths. More recent work has included making extensions to the HM approach to determine 
the linewidths when bounded contributions to the timing and phase jitter can be ignored ||5l, 
developing more elaborate models of pulse propagation in the laser Em, including technical 
noise contributions |I8][12|, and including the effects of gain dynamics ifTOl . A major unknown 
in all theoretical efforts to date was the strength of couplings between pulse parameters. The- 
oretical predictions carry large uncertainties and depend on the details of the laser design, in- 
cluding dispersion management. Recently, we quantitatively measured the linear response of 
the pulse energy, the central frequency, the round-trip gain ifTSl . the timing and the phase lfT4l 
of a mode-locked Ti:sapphire laser The motivation for that work was to quantitatively predict 
the SE-limited noise properties of the laser, the subject of this paper. 

When considering the shape of the comb lines of a mode-locked laser, there is a fundamental 
distinction between the case of a free-running laser and the case of a laser locked to an exter- 
nal oscillator In the former case, when the noise source is SE noise or any other white noise 
source, the central time and the phase of the mode-locked pulse undergo a random walk, and the 
comb lines have a stationary shape. In the latter case, again assuming white noise, the central 
pulse time and the phase are bounded, and there is no stationary line shape; the measured line 
width is inversely proportional to the measurement time. However, the phase noise spectrum 
of each comb line is stationary, and that spectrum — not the frequency spectrum — is physi- 
cally meaningful. Any clock or frequency measurement system consists of an oscillator and a 
counter ifTSl . The clock's performance depends on the frequency noise of the oscillator and the 



phase noise of the counter. Virtually all theoretical calculations to date of the noise properties of 
passively mode-locked lasers have focused on the frequency spectrum of a free-running laser, 
but in modern time and frequency metrology applications, these lasers are part of the count- 
ing system [4|, so it is their phase noise spectrum after they are locked to an oscillator that is 
important. Here, we calculate the SE-limited noise properties of a Ti:sapphire laser, comparing 
the comb line shape of the free-running laser and the residual phase noise of comb lines of the 
phase-locked laser, using experimentally-derived parameters for the linear response of the laser 
to noise perturbations. 

2. Linear response and noise drivers 

The laser is characterized by the pulse energy w ~ Weq + Aw, the gain g = ^eq + ^g, the central 
frequency flj ~ flJeq + ACT, the central pulse time T = Teq + At, and the phase 9 = 0eq + A0, 
where Xeq denotes the equilibrium value of each quantity and Ax the fluctuation about equilib- 
rium. Defining Tr as the round trip time, we have the repetition rate /rep = 1/7r. Without loss 
of generality, we may choose Teq = and 0eq = 27r/o7R, where /o is the equilibrium carrier- 
envelope offset frequency. We define the phase with respect to a common reference rather 
than the pulse envelope. Consequently, is related to the carrier-envelope phase slip AOcs by 
A0 = — A0ce + CTAt. The linear response of the laser to perturbations is JS] 

dv 

_A.v + S, (1) 

dT 

where v = (Ag,Avv',Anj,AT,A0)' (f denotes the transpose), S is a vector of noise sources for 
each parameter, and A is a matrix of coefficients that describe the linear response of each pa- 
rameter to changes in itself or the others. Once A and S are known, one can calculate the timing 
and phase jitter 13] , which leads to the optical frequency comb line shapes |f5| . The noise sources 
perturb all five parameters, and also add energy into a continuum of dispersive modes, but these 
modes must be damped in order for a passively mode-locked laser to be stable. Consequently, 
the noise contributions to those modes do not affect the linewidths. 

To calculate the frequency or phase noise spectrum, we use the formalism of stochastic dif- 
ferential equations lfT6l . The second-order moment equations for the laser parameters are 

ii(Av;Av;) w^, 

^ ,' " = - E(^'*- +^Ji) i^^'k^^i) + Dij, (2) 

"^^ k.i 

where v; is one of the dynamical variables and {Si{T)Sj{T')) = Dij5{T — T') relates S to a 
matrix D of noise drivers. For the 5 dynamical parameters introduced above, there are thus 
15 differential equations describing the response of the second-order moments to the noise 
sources. The quantum noise drivers Djj were found using a generalized perturbation theory for 
Gaussian pulses in a dispersion managed soli ton laser |fT3l . Defining the cavity lifetime Tph, 
the fluorescence lifetime Tf , the equilibrium photon number A'ph.eq, the equilibrium upper state 
population Af2,eq, and the pulse width fp, they are 

2 1 ^^^eq „ „ A^ph,eq 1 /iCJeq 

D,„, =2Weq ;;- — , D„g = Dg„. =Weqgeq 
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where F{tp) = C\{tp)-C,{t,)C2{tp) + (3/4)C2(fp), Ci(fp) = (1 +3.A,1)/(1 +5s/t% C2{t„) = 
(1 +i/fp)/(l +5i/d), and 5 is a dispersion management parameter. The resuhs above are sim- 
ilar to previously used quantities derived from ordinary soliton perturbation theory fT\. While 
most of the drivers differ by small factors, the energy and gain drivers are identical. 

For the laser studied in [HI and HU, Weq = 55 nJ, §eq = 0.054, CJeq = 2.3 x 10^^ rad/s, 
Tph = 0.1 /is, A'2,eq/A^ph,eq ~ 13, and s/tp = 0.05, which leads to F{tp) = 0.60. We use the 
experimentally-measured values from ITSJI and lfT4l for the elements of A. One could solve 
Eq. ^ in the time domain numerically, but there are practical difficulties arising from the wide 
range of time scales involved in the dynamics. We instead found the moments in the frequency 
domain, using the formal solution 

(Ay,„Av;)(») =(-/«! - A)-' • D • [(/col - A)"']', (4) 

where I is the identity matrix. We then used the inverse Fourier transform to bring this solution 
into the time domain symbolically. 

3. Free-running lineshape 

Taking into account only the effects of timing and phase jitter (amplitude and central frequency 
noise both contribute to the Unewidth, but are in practice negligible, c.f. ||2l), the power spectral 
density of the «* comb line is 

/„(©) - j|"fl'rcos(a)r)exp |-iii2([AT(r)]2) +n„(AT(r)A0(r)) - i([A0(r)]2)| , (5) 

where o is the Fourier frequency with respect to the center of the comb line O),,, and ii„ = 
O,, — CTeq. See Fig. [T]for plots of the free-running noise properties. With no feedback on /^ep 
and /o, the second-order moments involving only t and do not evolve to equilibrium values at 
large T , but rather grow with time. Each can be expressed as yT + P (T), where 7 is a constant 
and j3(r) -^ as T ^ °o. The calculated single-sided shape for two comb lines, one in the 
center and one on the wing of the spectrum, is shown in Fig. [T^. The dynamics of the noise on 
short time scales arising from the residual oscillatory part j3(r) of the second-order moments 
lead to features on the wings of the spectrum, as shown in Fig. [T^. The most notable feature 
is the peak at 400 kHz, the relaxation oscillation frequency. The width of this peak depends 
on the strength of the effective saturable absorber that mode-locks the laser, which is strongly 
influenced by the pump power IIT3JI and can become overdamped under certain conditions. 

The line is to a very good approximation Lorentzian. The FWHM linewidth of the «* comb 
line is simply jn = Yee —2Q.„YeT + ^nYTT, where jxy- is the asymptotic slope of each correspond- 
ing term in the exponent of Eq. Q. The linewidth is plotted as a function of i2„ in Fig. [T]?. The 
noise in the phase A0 is dominated by energy fluctuations coupling into phase noise through the 
Kerr effect. The Schawlow-Townes linewidth, the linewidth if the phase noise is driven by Dgg 
alone, is more than two orders of magnitude smaller. The noise in the time At is dominated 
by fluctuations in the central frequency coupling into timing jitter through cavity dispersion. 
These central frequency fluctuations are driven mostly by the direct term Dcjaj. The noise is 
dominated by timing jitter over most of the spectrum, suggesting the importance of minimiz- 
ing the cavity dispersion for comb applications. Also note that the correlation term (A0At) 
leads to a shift in the frequency of minimum linewidth from the central frequency. This shift is 
especially important for the Kerr effect, which affects both phase and timing IfTTl . 

4. Feedback, the phase noise spectrum, and frequency uncertainty 

In metrological applications of femtosecond frequency combs, the comb is phase-locked to a 
frequency standard, and this has a profound effect on the lineshape. The moments involving 
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Fig. 1. The calculated free-running comb lineshape. (a) Power spectrum I„{co) for comb 
lines near the central frequency (black) and 200 THz from the central frequency (red), (b) 
Linewidth y as a function of Q.„ , the offset from CJeq . The red curve is the result with the full 
set of A coefficients measured in 1131 and 1141 . The black curve is the result when timing 
jitter caused by dispersion is neglected (by setting A^cj = 0). The blue curve is the result 
when, in addition, timing and phase jitter driven by intensity fluctuations is neglected (by 
setting Atv,. = and Ag^. = 0). 



only T and 6 are now damped, and their time domain solutions now take the form a + 77(r), 
where a is a constant and ri{T) ^ as T ^ °o. The first term leads to a delta function in 
frequency with an area proportional to exp(— a). This part of the comb line is coherent, which 
makes it useful in metrology applications. The second term coiTesponds to residual phase noise, 
the shape of which depends on the laser dynamics and the details of the feedback system. In 
metrology applications, the residual phase noise leads to uncertainty in frequency measure- 
ments ifTSl . 

In optical frequency metrology H or optical clock development and comparison llT8l[T9l . 
/o is measured and locked using self-referencing and /rep is locked to a microwave frequency 
standard or to a heterodyne beat between a comb line and a laser line locked to an optical 
atomic transition. Because it depends on the phase with respect to the envelope, /o depends 
on both T and 6. In general, the phase of a comb line depends on both parameters as well, 
so each feedback system, for both /o and /rep, will measure and act on fluctuations in both T 
and 9. If, for simplicity, one assumes that the feedback mechanism — a cavity miiTor on a 
piezoelectric transducer or an acousto-optic modulator controlling the pump power — acts on 
A/o and A/„, rather than some other linear combination of At and A0, one can proceed by 
defining a linear transformation which goes from the {t, 9} basis to the {/o,/i} basis. Each 
feedback loop acts on the latter pair of variables, and the inverse transformation determines how 
the feedback affects T and 9. We have modeled a feedback system with parallel proportional 
and integrator stages because it is simple enough that analytical solutions can still be found. The 
parameters we choose give a servo bump at about 1 kHz. We assume a simple radio frequency 
/rep lock rather than a lock using an optical heterodyne beat. We have found that shot noise in 
the detection of beats used for locking the laser introduces a negligible amount of noise through 
the feedback loops because the bandwidth of a typical feedback system is approximately two 
orders of magnitude smaller than /rep. We note that when the feedback bandwidth becomes 
comparable to the dynamical frequencies in the laser, the laser dynamics participate in the 
feedback mechanism ifTOl . That is not the case in the system that we considered. 
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Fig. 2. Calculated phase noise properties of the locked comb. Phase noise spectrum for 
comb lines near C7eq (black), 100 THz away (red), and 200 THz away (blue). 



For the «* comb line, the phase noise spectrum is 

5,t(a)) = iii2([AT(to)]2)-n„(AT(to)A0(to)) + ^([A0(to)]2). 



(6) 



We show this solution for several comb lines in Fig. |2]for the given laser and feedback parame- 
ters. While shot noise does not add much noise to the comb lines through the feedback, it does in 
general contribute strongly to the phase noise of a heterodyne beat, but we do not consider that 
here. The phase noise spectrum contains all the information needed to determine the uncertainty 
in a frequency measurement. There has been recent work on how best to estimate uncertainty 
from measurements using modern frequency counters, which use interpolation to increase their 
accuracy over older style counters 12011211 . Counting of zero-crossings is well-described by 
the Allan variance, whereas the uncertainty in measurements using modern counters is better 
described by the triangle variance 1211 . We have calculated both uncertainty estimates using 
the SE-limited phase noise spectrum and found that, in the center of the spectrum where the 
noise is smallest, the Allan (triangle) deviation is 8x10^'^ (8x10^^') for 1 second averaging 
time. The large difference between the two arises from the triangle variance's stronger rejection 
of high frequency (compared to the gate time) noise. For a comb line 200 THz away from the 
center of the spectrum, we find that the Allan (triangle) deviation is 1 x 10^^^ (1 x 10^^") for 1 
second averaging time. Note that these values depend strongly on the feedback parameters. 

5. Conclusion 

Combs have already been shown to be comparable with 1-s Allan deviation on the order of 
10^'^ fractional uncertainty ll22l . The calculations presented here are of the same order of 
magnitude, but we caution that our calculation is specific to the laser studied in ifTSI and lfT4l . 
using the simple feedback model we considered. We note that our approach applies to technical 
noise sources as well, if one adds them to the driving terms contained in D. Using that informa- 
tion, coupled with experimentally-derived parameters for A and a detailed and accurate model 
of the feedback system, it should be possible to make quantitative predictions of uncertainty 
in frequency measurements. A program focused on optimizing the linear response of the laser 
system could be used to produce better laser designs for optical frequency metrology and other 
applications of femtosecond frequency combs. 
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